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ABSTRACT 


This  paper  presents  a simple  approach  for  the  study  of 
the  following  two  types  of  multiserver  finite  queueing  systems 
each  with  two  classes  of  customers  demanding  service  with  different 
arrival  and  service  rates:  (1)  Customers  belonging  to  one  class 
have  priority  over  the  other  and  the  number  of  low  priority 
customers  that  may  be  present  in  the  system  at  any  time  is  restricted. 
Customers  blocked  from  entering  service  at  any  time  are  cleared  from 
the  system.  (2)  Customers  belonging  to  one  class  are  cleared  when 
blocked  (loss  system)  and  customers  belonging  to  the  other  class  are 
buffered  (delay  system)  under  such  circumstances.  Furthermore,  the 
number  of  customers  belonging  to  the  loss  system  demanding  service 
at  any  time  is  restricted. 


I 
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INTRODUCTION 


Consider  the  following  two  classes  of  multi-server  finite 
queueing  systems: 

1)  Multi-server  loss  systems  with  two  classes  of  customers 
with  different  arrival  and  service  rates.  Customers  belong- 
ing to  one  class  have  priority  over  the  other  and  the  number  of  low 
priority  customers  that  may  be  present  in  the  system  at  any  time  is 
restricted  by  a number  K (K  < s,  number  of  servers).  These  systems 
will  be  denoted  as  A/B^.B^/s/s. 

2)  Multi-server  systems  with  two  classes  of  customers  in  which 
customers  belonging  to  one  class  are  cleared  when  blocked  (loss  system) , 
and  customers  belonging  to  the  second  class  are  buffered  under  such 
circumstances  (delay  system) . The  two  classes  have  different  arrival  and 
service  rates.  As  before,  tne  number  of  customers  belonging  to  the 

loss  system  that  may  demand  service  at  any  time  is  restricted  by  K. 

We  denote  this  class  of  systems  by  A/B, ,B„/s/N. 

Such  queues  are  common  in  the  study  of  telecommunications  networks  in 
which  several  classes  (such  as  voice  and  data)  share  access  and  traffic 
line  facilities.  Under  purely  Markovian  assumptions,  the  simplest  of 
such  systems  have  been  analyzed  by  Cohen  [2].  However,  when  the 
number  of  servers  is  greater  than  one,  the  magnitude  of  the  problem 
increases  considerably  as  illustrated  in  Kotiah  and  Slater  [3] , where 
a two  server  system  of  this  type  has  been  analyzed.  For  larger  values 
of  s,  even  though  an  exact  analysis  is  possible  in  principle,  in 
practice  there  are  severe  computational  limitations  (as  also  pointed 
out  by  Bhat  and  Fischer  [1]).  An  approximate  technique  to  determine 
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some  of  the  system  performance  measures  has  been  presented  in  Bhat 
and  Fischer. 

Here  we  are  concerned  with  the  steady  state  marginal  distributions 
of  each  class  from  which  one  could  easily  deduce  the  expected  values 
and  variances  of  number  of  customers  and  the  respective  blocking  pro- 
babilities. Bearing  in  mind  the  utility  of  the  derived  results  and 
keeping  the  generality  of  the  above  description  to  a large  extent, 
we  shall  assume  the  following  common  characteristics  for  the  two 
classes  of  multi-server  systems  mentioned  above: 

(1)  The  two  classes  of  customers  arrive  in  a single  recurrent 
process  with  mean  (1/X).  An  arriving  customer  belongs  to  one  of  the 
two  classes  with  constant  probabilities  p and  1-p  respectively. 

(2)  Both  classes  of  customers  have  exponential  service  times 
with  rates  and  respectively. 

(3)  There  are  s servers  in  parallel  and  if  buffering  is 
allowed,  no  more  than  n customers  are  allowed  in  the  system  at  one 
time. 

(4)  An  arriving  class  2 (low  priority)  customer  can  not  enter 
the  system  if  K of  them  are  receiving  service  already  and/or  the 
buffer  has  some  clas^  1 customers  waiting. 

With  these  simplifications,  the  two  classes  of  queueing  systems 
may  now  be  represented  as  GI/M^,M2/s/s  and  GI/Mj^,M2/s/N  respectively. 

In  the  above  description  of  the  problem,  the  queue-length  process, 
when  restricted  to  the  arrival  epochs,  describes  a discrete  time 
Markov  chain.  Clearly,  the  measures  of  interest  can  be  easily  deduced, 
if  we  obtain  the  stationary  distrlbutioh  of  this  imbedded  Markov  chain. 
However,  in  reality,  for  moderate  to  large  systems,  the  number  of  states 
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[ (K+1) (N+l-K/2)  for  GI/M^.M^/s/N,  N > s]  is  going  to  be  too  large 
to  make  an  exact  numerical  approach  for  the  determination  of  the 
stationary  distribution  of  the  Markov  chain  possible.  Furthermore, 
because  of  the  cumbersome  nature  of  the  expressions  involving 
binomial  coefficients  for  the  one  step  transition  probabilities  of 
the  Markov  chain,  generation  of  the  transition  probability  matrix 
(TPM)  itself  might  pose  considerable  numerical  problems.  This  will 
be  exemplified  in  the  following  section.  Therefore,  alternate 
techniques  are  needed  to  provide  at  least  approximate  values  for 
system  performance  measures. 

In  this  paper  we  develop  approximation  techniques  for  determining 
the  first  two  moments  of  system  sizes  and  the  probabilities  of  blocking 
for  each  of  the  customer  classes.  Exact  analysis  is  used  to  ascertain 
the  adequacy  of  our  heuristic  procedure.  Based  on  our  computational 
experience,  the  approximation  provides  better  results  for  the  system 
GI/M^,M2/s/s  than  for  the  system  GI/Mj^,M2/s/N.  Nevertheless,  in  the 
absence  of  an  analytical  justification  these  conclusions  are  restricted 
to  the  class  of  systems  used  in  the  analysis.  Incidentally,  since  the 
expressions  for  the  one  step  transition  probabilities  in  the  case  of 
GI/Mj^,M2/s/N  are  very  complicated  for  ntmierical  work,  a justification 
for  the  procedure  should  come  from  an  elaborate  simulation  study  of 
these  systems,  which  in  itself  is  a major  problem. 

In  the  systems  of  interest,  if  * ^2  * numerical  methods 
developed  in  Raju  [4]  and  Raju  and  Bhat  [5]  would  directly  yield  the 
measures  of  interest  since  there  is  no  need  to  distinguish  the 
customers  when  once  they  join  the  system.  In  an  indirect  way,  the 
procedure  used  in  this  paper  heavily  utilizes  this  situation.  For 
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further  development  we  introduce  the  following  notations: 

A(x)  : Inter-arrival  time  distribution  (0  < x <«>) 

: Traffic  intensity  of  class  i (i  = 1,  2)  customers 

( = for  i=l  and  = Xd-p)/^!^  for  i“2) . 

i til 

: Number  of  class  i customers  at  the  n'"  arrival  epoch 

in  the  system. 

til 

Q : Number  of  customers  in  the  system  at  the  n arrival 
n 

epoch  ( = + Q^) . 

n n 

Define, 

■ ' ‘C  - Cl  • 'loi  " *• 


Thus  represents  the  transition  probability  of  the  process 

1 2 

(Q  »Q  ) from  i-^r  and  i-*-t  respectively  during  an  interarrival  time 
n n 

period.  The  TPM  of  the  imbedded  Markov  chain  has  a's  for  its  elements. 
With  this  background,  we  shall  next  consider  the  approximations  for  the 
respective  systems. 


APPROXIMATION  FOR  GI/M^.M^/s/s  SYSTEM: 

Recall  that  an  arriving  customer  is  of  class  1 with  probability  p 
and  of  class  2 with  probability  (1-p) . For  the  sake  of  simplicity,  we 
let 

£(.)  • 1 - 
8(n)  - 1 - *-"''2* 

and 

h(m,n)  - e ^”'^1  (0  < x < ») 

We  have  for  this  system 
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(^)  r°  f(i+l-r)  g(j-t)  h(r,t)  dA(x) 
ij  r t Q 

+ (l-p)  (^)  / f(i-r)  g(j+l-t)  h(r,t)  dA(x) 

^ 0 

(0<i  + j<s,  0'<j<K,  0<r<i  + l,  0<t<j+l) 

“c,  •«  /**  f(s-j-r)  g(j-t)  h(r,t)  dA(x) 

®**J  ♦ J 9^^  ^ ^ 0 

(j  < K,  r < s - j,  t < j) 

and 

(^)  r f(i+l-r)  g(K-t)  h(r,t)  dA(x) 

Xlx^lTu  X L Q 

+ (l-p)  (^)  (^)  P f(i-r)  g(K-t)  h(r,t)  dA(x) 

(i  + j < s,  r < i + 1.  t < K) 

The  cumbersome  nature  of  the  expressions  is  evident  and  moreover,  it 
does  not  seem  possible  to  get  around  the  problem  of  direct  evaluation 
of  binomial  coefficients.  Thus,  even  for  moderate  values  of  s and/or  K, 
the  generation  of  TPM  might  turn  out  to  be  expensive  and  unreliable. 

However,  if  is  clear  that  {Q^  : n - 1,  2,  3,  . . . } is  a 

Markov  chain  imbedded  in  the  queue-length  process.  Additionally  assuming 
K “ 3 in  the  above  system,  we  obtain  a system  of  the  type  GI/M/s/s.  Let 

Pfj  - P - j I - 1}  (0  < i,  j < s). 

we  have 

p r dA(x) 

J 0 

+ (l-p)  (J)  jf*  (l-e-’^")^-J  e-^^"  dA(x) 

(i  < s,  j < i + 1) 

(!)  r dA(x)  (j  < s,  i - s) 

J 0 
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The  resulting  matrix  P = 1|  P^, ^ ||  is  almost  left  triangular  (or  lower 
Hessenberg)  and  can  be  easily  generated.  Denote 


C.J  - dA(x) 


(0  < j < s) 


After  some  algebraic  manipulation,  we  can  establish  that  P can  be 
generated  in  two  stages  through  the  following  recursive  relations. 


First  Stage 


P 


i,i+l 


pa 


i+1 


and 


Pii  * pw-p”  '“i  - Vi’ 


(0  < i < s) 


i+1  1+1 

i-j+1  ^i-l,j  “ i-j+1  ^ij+1 


(1  < i < s,  j < i) 


P 


sj 


(0  < j < s) 


Second  Stage 

Starting  from  i^s-l,  s-2,...,l  and  for  j < i 


’ = p + ^ P p 

ij  ^ij  p ^1-lj 


Then  reset  the  diagonal  entries  to 


+ (1-p)  ot^  (0  < 1 < s) 

Remark:  It  nay  be  noted  that  if  we  set  p=l  in  1 , we  get  the  expres- 

sions for  the  ordinary  GI/M/s/s  system.  Using  the  first  stage  of  the 
above  recurrence  relations  with  p*l,  we  can  very  efficiently  generate 
the  TPM  for  this  system. 

The  objective  is  to  approximate  the  measures  of  interest  for  the 


7 


system  GI/M^,M^/s/s  by  the  corresponding  measures  of  two  different 
systems  of  the  type  GI/M/s/s  for  which  the  characteristics  can  be 
easily  determined  by  the  methods  developed  in  Raju  [4].  For  the 
sake  of  clarity,  we  shall  denote  this  approximating  system  as  GI*/M/s/s. 
Given  below  is  the  approximation  scheme  for  GI/M^jM^/s/s  system. 

The  Heuristic: 


Initialize: 

N, 


[K  • P2/(Pi  + 


= s - 


Step  1: 


Solve  GI*/M^/N^/N^  (i.e.,  in  ( 1 ) , put  y s =• 

obtain  N*  = E[Busy  Servers]  in  the  above  system 


Set 


[s  - N^l 


Restart:  = mlnCN^.K) 


Step  2: 

Solve  GI*/M2/N2/N2  (in  (1)  , put  y * y^ , s ” and  inter- 
change p and  1-p) 

Obtain  the  stationary  distribution  {1T2  “ 1^2^  : 0 < j < N2}  of 
the  imbedded  Markov  chain.  Let  E2  and  V2  be  the  expected  value 
and  variance  of  the  distribution  respectively. 

Set  » (s-E2)'‘‘ 

Step  3: 

Solve  GI*/Mj^/N^/N^  (in  (1)  , put  y - y^  , s - N^) 

Obtain  the  stationary  distribution  : 0 < j < 

of  the  imbedded  Markov  chain.  Let  and  be  the  expected 
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value  and  variance  of  the  distribution  ^ respectively. 

Step  4: 

Stop  (or  go  to  Restart  with  N2  = (s-E^^)^). 

The  measures  of  interest  for  the  system  GI/M^.M^/s/s  are  now  approximated 

by  E^,  and  which  represent  for  class  i (i=l,2),  the  expected  num- 

i 

ber  of  customers,  variance  of  the  number  of  customers  and  probability  of 
blocking  respectively. 

Remark : If  one  wishes  to  improve  upon  the  results  obtained  after 
the  first  iteration,  one  could  go  to  Restart  from  step  4.  However,  we 
do  not  have  enough  evidence  to  suggest  that  such  an  attempt  would  always 
yield  better  results. 

APPROXIMATION  FOR  GI/M^.M^/s/N  SYSTEM 

The  basic  idea  behind  .^our  procedure  is  to  approximate  the  characteristics 
of  interest  for  the  system  GI/M^,M2/s/N  by  the  corresponding  measures  of 
two  different  systems  of  the  type  GI/M/s/N  for  which  the  measures  can  be 
readily  obtained  by  the  methods  developed  in  Raju  [4].  For  the  sake  of 
clarity  we  shall  denote  the  approximating  system  as  GI  /M/s/N. 

The  Heuristic; 

Initialize: 

S2  “ [K  • , N2  = K 

Step  1: 

Solve  GI*/M2/s2/N2, 

it 

obtain:  Sj^  = E [number  of  customers  in  the  system] 

Step  2: 

Solve  GI  /M^/sj^/Nj^,  where 
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= " - Sr  N,  = N - . 


Obtain  the  stationary  distribution  0 < j < of  the 
imbedded  Markov  chain,  hence  get  E^,  V^,  the  expected  value 
and  variance  of  the  stationary  distribution  respectively. 


Obtain  s^  = E [busy  servers]  due  to  the  above  system. 


Step  3: 


Solve  GI  /M2/S2/N2,  where 

* + 

= [min  [K,  s - s^]  ] . 

Obtain  the  stationary  distribution  ’ 0 ^ j ^2)  of  the 
imbedded  Markov  chain.  Let  E^,  ^2  expected  value  and 

variance  of  this  distribution  respectively. 


Step  4:  STOP. 

The  measures  of  interest  for  the  system  GI/M^^.M^/s/N  are  now  approxi- 
mated by  E^,  and  which  represent  for  class  i(i=l,  2)  the  expected 

number  of  customers,  variance  of  the  number  of  customers  and  probability 
of  blocking  respectively. 


NUMERICAL  ILLUSTRATIONS: 

A.  The  System  GI/M^,M2/s/s 

For  the  purpose  of  illustrating  the  approximation  developed  for  the 
system  GI/M^,M2/s/s  in  the  precedding  section,  we  present  below 
tables  containing  exact  measures  of  performance  for  the  system  and  their 
respective  approximations  obtained  through  our  heuristic  under  special 
cases.  It  is  assumed  that  the  time  intervals  between  successive  arrival 
epochs  have  an  Erlangian  distribution  (Gamma  distribution)  with  a prob- 
ability density  function 


f(x) 


-mAx 

e 


(mX)*"  x”' 
(m-1) ! 


1 

— dx 


m*l , 2 , 3 , 


A > 0 . 
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The  arrival  rate  is  given  by  X.  The  following  notations  are  used. 

System  parameters: 

s : system  capacity  and  the  number  of  servers. 

K : system  capacity  for  class  2 access 
p ; probability  of  class  1 arrival 

X : arrival  rate 

m : number  of  phases  in  the  Erlangian  distribution 
c^  : coefficient  of  variation  of  arrival  process 

: 1.0/t^ 

; service  rate  for  class  1 customers 

y2  • service  rate  for  class  2 customers 

Measures  of  performance:  E(Q^)  : mean  number  of  class  1(=1  or  2) 

customers  at  arrival  epochs. 

V(Q^)  : variance  of  number  of  class  i(=l  or 
2)  customers  in  system  at  arrival 
epochs 

PB^  : probability  of  blocking  for  class  i 
(i=l,  2). 

Approximating  systems  parameters: 

N^:  system  capacity  and  the  number  of  servers 

for  class  i (=1,  2) 

X^:  arrival  rate  for  class  i (»1,  2) 

« pX  if  i=l,  ” (l-p)X  if  i=2. 

The  interarrival  times  for  class  i («1,  2)  system  are  assumed  to  have 
an  Erlangian  distribution  with  probability  density  function 


f^(x) 


-m  X.x 
e i 


...  m m™  1 
(mX^)  X 

(m-1) ! 


> 0 


dx 
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The  measures  of  performance  of  the  system  GI/M^,M2/s/s  are  approxi- 

* ★ 

mated  by  the  measures  of  performance  of  GI  /M^/N^/N^  and  GI 
systems.  In  Tables  la  and  lb,  for  each  set  of  variables  the  exact 
measures  for  the  system  GI/M^^.M^/s/s  are  presented  and  the  corresponding 
approximations  obtained  from  the  approximating  systems  appear  immediately 
below.  For  instance,  the  digits  inside  the  braces  following  the  number 
of  servers  in  the  system  GI/Mj^,M2/s/s  represent  and  N2  respectively. 

Based  on  our  numerical  work  the  following  observations  can  be 

made. 

1.  The  approximation  performs  well  under  low  traffic  conditions 
which  are  combinations  of  the  arrival  rate,  service  rate  and 
system  capacity. 

2.  The  effect  of  the  approximation  is  more  pronounced  on  the 
variance  of  the  number  of  customers  in  the  system. 

3.  The  effect  of  the  approximation  is  more  pronounced  under 
higher  variability  conditions  of  the  arrival  process. 

B.  The  System  GI/Mj^,M2/s/N 

Because  of  the  analytically  intractable  nature  of  the  one-step 
transition  probability  expressions  (obtained  by  extending  the  expres- 
sions given  in  (1)  for  cases  s < i,j  < N) , the  exact  analysis  of  this 
system  is  extremely  cumbersome  and  time  consuming.  However,  if  we 
assume  exponential  inter-arrival  times  for  each  of  the  two  classes 
of  customers,  the  desired  characteristics  of  the  above  system  can  be 
obtained  by  solving  the  balance  of  state  equations.  (See,  Bhat  and 
Fischer  [1]).  Therefore,  numerical  comparisons  for  this  system 
approximating  procedure  shall  be  restricted  to  Poisson  arrivals. 
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s / 

ni» 

M2  / 

S / S SYSTEM 

s 

R 

L*i 

U1 

U2 

CQV 

E(OII 

E(021 

V(011 

VIC2) 

EBl 

M82 

12 

4 

6.000 

2,000 

.400 

.354 

I.Q4 

3.49 

.89 

• 51 

•0000 

• 3607 

^1 

- 

1.04 

3.49 

.89 

• 51 

• 0000 

• 3607 

12 

4 

6.000 

2.000 

*400 

.447 

1.05 

3.49 

• 91 

.52 

.0000 

• 3614 

9, 

• 

1.05 

3.49 

• 91 

• 52 

• 0000 

• 3614 

12 

4 

6.000 

2.000 

.400 

.707 

1.11 

3.48 

1.02 

.54 

• 0000 

• 3641 

- 

1.11 

3.48 

1.02 

• .54 

• 0000 

• 3641 

12 

4 

6.000 

2.000 

.400 

1.000 

1 .20 

3.48 

1.20 

• 58 

• 0000 

• 3663 

9, 

- 

1.20 

3.46 

1.20 

.58 

.0000 

• 3683 

12 

4 

6.000 

2.000 

• 400 

1.545 

1.38 

3.47 

1.49 

• 65 

• 0000 

.3778 

9» 

“ 

1.38 

3.47 

1.49 

• 65 

• 0000 

.3778 

12 

4 

6.000 

2.000 

• 400 

2.047 

1.54 

3.48 

1.73 

• 69 

• 0001 

• 3875 

9# 

1.54 

3.48 

1.73 

• 69 

• 0000 

.3875 

12 

4 

12. 000 

2.000 

• 400 

• 354 

2.23 

3.76 

1.87 

.24  . 

• 0001 

.4728 

9, 

* 

• 2.23 

3.76 

1.67 

fi** 

• 0000 

.4728 

12 

4 

12.000 

2.000 

.400 

.447 

2.24 

3.76 

1.92 

• 25 

• 0002 

.4730 

9, 

- 

2.25 

3.76 

1.92 

• 25 

.0000 

.4730 

12 

4 

12.000 

2.000 

• 400 

.707 

2.30 

3.76 

2.09 

• 25 

• 0003 

.4737 

9, 

41 

- 

2.30 

3.76 

2.10 

• 25 

• 0001 

.4737 

12 

4 

12.000 

2.000 

.400 

1.000 

2.40 

3.75 

2.37 

.27 

• 0008 

.4748 

9, 

41 

" 

2.40 

3.76 

2.39 

.27 

• 0003 

• 4748 

12 

4 

12.000 

2.000 

.400 

1.545 

2.60 

3.75 

2.86 

.29 

• 0022 

.4777 

9f 

4) 

2.61 

3.75 

2.94 

.29 

• 0008 

.4777 

12 

4 

12.000 

2.J00 

.400 

2.047 

2.8? 

3.75 

3.32 

• 32 

• 0040 

.4810 

9, 

41 

2.83 

3.75 

3.43 

• 32 

.0016 

.4809 

12 

4 

18.000 

2.0CO 

• 400 

.354 

3.40 

3.84 

2.73 

• 16 

• 0030 

.5138 

9, 

4) 

- 

3.42 

3.84 

2.81 

• 16 

• 0010 

.5137 

12 

4 

18.000 

2.000 

.400 

.447 

3.42 

3.84 

2.78 

• 16 

• 0033 

.5138 

9» 

41 

- 

3.43 

3.84 

2.87 

• 16 

.0012 

.5138 

12 

4 

18.000 

2.000 

.400 

.707 

3.46 

3.84 

2.97 

• 16 

• 0046 

.5141 

9# 

41 

• 

3.49 

3.84 

3.08 

• 16 

.0018 

• 5141 

12 

4 

18.000 

2.000 

.400 

1.000 

3.54 

3.84 

3.25 

.17 

.0070 

.5147 

9* 

41 

- 

3.57 

3.84 

3.42 

.17 

• 0031 

• 5146 

12 

4 

18.000 

2.000 

.400 

1.545 

3.72 

3.84 

. 3.79 

.18 

• 0124 

• 5160 

9, 

41 

- 

3.77 

3.84 

4.08 

• 19 

• 0062 

• 5159 

12 

4 

NlS.OOO 

2.000 

• 400 

2.047 

3.91 

3.83 

4.24 

.20 

•0183 

.5176 

9* 

41 

3.99 

3.84 

4.66 

• 20 

• 0100 

.5175 

13 


s / 

flit 

M2  / S 

/ S SYSTEM 

K 

t*M 

Cl 

C2 

cov 

EOU 

ff  02) 

V(Q2  ) 

V(Q2) 

Pff2 

2 

3.000 

l.CCO 

6.290 

.394 

1.03 

• 11 

• 68 

.10 

• 0007 

.0024 

)f 

2) 

” 

1.03 

• 11 

• 86 

• 10 

• 0009 

• 0014 

2 

3.000 

l.CCO 

6.290 

.447 

1.09 

• 13 

• 90 

• 12 

• 0009 

.0034 

5# 

Zl 

• 

1.09 

• 13 

• 91 

• 12 

• 0006 

• 0022 

2 

3.000 

l.OCO 

6.290 

.707 

1.10 

• 19 

1.00 

.17 

• 0021 

.0096 

S> 

2) 

- 

1.10 

.19 

1.01 

• 18 

.0012 

.0069 

2 

3.000 

l.OCO 

6.290 

1.000 

1.16 

• 28 

1.14 

• 26 

• 0092 

• 0244 

21 

- 

1.19 

• 26 

1.16 

• 26 

• 0029 

• 0197 

2 

3.000 

1.000 

6.290 

1.549 

1.34 

• 39 

1.34 

• 32 

.0106 

.0412 

S» 

2) 

• 

1.36 

• 36 

1.40 

.33 

• 0052 

• 0304 

2 

3.000 

1.000 

6.290 

2.047 

1.48 

• 40 

1.49 

.36 

• 0158 

• 0538 

2) 

1.51 

• 41 

1.98 

• 37 

• 0081 

• 0361 

2 

6.000 

l.OOO 

6.290 

.394 

?.ll 

.33 

1.53 

• 28 

.0218 

.0472 

5t 

2) 

• 

2.15 

• 39 

1.59 

• 29 

• 0151 

.0191 

2 

6.000 

1.000 

6.290 

.447 

2.11 

.35 

1.94 

.29 

• 0236 

.0918 

S» 

21 

• 

2.16 

• 37 

1.62 

• 31 

• 0160 

• 0223 

2 

6.000 

1.000 

6.290 

.707 

2.13 

• 40 

1.59 

• 34 

• 0309 

.0697 

5> 

21 

• 

2.19 

• 43 

1.71 

• 36 

.0194 

• 0354 

2 

6.000 

1.000 

6.290 

1.000 

2.19 

• 48 

1.69 

• 41 

• 0414 

.0972 

5# 

2) 

• 

2.29 

.52 

1.64 

.44 

• 0290 

.0571 

2 

6.000 

1,000 

6.290 

1.5*5 

2,24 

.57 

1.73 

• 46 

• 0573 

• 1311 

9# 

21 

• 

2.38 

• 63 

2.03 

• 50 

.0359 

• 0814 

2 

6.000 

1.000 

6.290 

2.047 

2.34 

• 62 

1.86 

• 49 

• 0705 

.1559 

9* 

21 

2.92 

• 70 

2.19 

.94 

.0455 

.0978 

2 

9.000 

1.000 

6.290 

.394 

2.92 

.49 

1.54 

• 38 

.0719 

.1343 

9> 

21 

• 

2.97 

• 96 

1.68 

• 43 

• 0930 

.0546 

2 

9.000 

1.000 

6.290 

,447 

2.81 

• 50 

1.54 

« .39 

.0736 

.1390 

9» 

21 

• 

2,98 

.59 

1.69 

• 44 

• 0941 

.0590 

2 

9.000 

1.000 

6.290 

.707 

2.79 

• 94 

1.97 

• 43 

• 0815 

.1961 

9f 

2) 

• 

2.99 

• 64 

1.75 

• 48 

.0989 

.0754 

2 

9.000 

1.000 

6.290 

1.000 

2,76 

.60 

1«60 

.47 

.0928 

• 1 804 

9* 

21 

• 

3.01 

.72 

1.89 

• 94 

• 0692 

• 1001 

2 

9.000 

1.000 

6.290 

1.549 

2.73 

• 66 

1.66 

.91 

• 1099 

• 2139 

9* 

21 

• 

3.08 

• 83 

2.00 

• 98 

• 0784 

• 1309 

2 

9.000 

1.000 

6.290 

2.047 

2.83 

.73 

1.74 

.93 

• 1241 

.2391 

9* 

21 

3.17 

• 90 

2.10 

• 60 

• 0907 

.1917 
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Tables  2a  and  2b  present  exact  measures  of  performance  and 
their  respective  approximations  obtained  through  our  heuristic 
procedure.  The  following  notations  are  used: 

Systems  Parameters 

N : System  capacity 
s : Number  of  servers 
K : System  capacity  for  class  2 access 
X : Arrival  rate  of  class  iC*!,  2)  customers 
: Service  rate  of  class  i(=l,  2)  customers 
Measures  of  Performance: 

E(Q^)  : Mean  number  of  class  iC^l,  2)  customers  in  the 
system. 

V(Q^)  : Variance  of  class  i(=l,  2)  customers  in  the 
system. 

PB^  : Probability  of  blocking  for  class  i(=“l,  2). 

For  each  set  of  system  parameters,  the  approximate  measures  of 
performance  obtained  through  our  heuristic  are  presented  right 
below  the  corresponding  exact  measures. 

A general  observation  that  can  be  made  is  that  the  approximation 
gives  poorer  results  under  heavier  traffic  conditions.  Also  for 
reasons  which  need  further  investigations,  the  performance  of  the 
heuristic  procedure  on  this  class  of  systems  is  poorer  than  that 
of  the  previous  one. 
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